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SUMMARY 


The present report is a study of the stability of laminar boundary- 
layer profiles on slightly curved walls relative to small disturbances, 
in the shape of vortices, whose axes are parallel to the principal direc- 
tion of flow. The result in an eigenvalue problem by which, for a given 
undistvirbed flow at a prescribed wall, the amplification or decay is com- 
puted for each Reynolds number and each vortex thickness. For neutral 
disturbances (amplification null) a critical Reynolds number is determined 
for each vortex distribution. The numerical calculation produces ampli- 
fied distixrbances on concave walls only. The variation of the dimension- 


less 


Uo^ 


— with respect to ad is only slightly dependent on the 
R 

shape of the boundary-layer profile. The numerical results yield informa- 
tion about stability limit, range of wave length of vortices that can be 
amplified, and about the most dangerous vortices with regard to the tran- 
sition from laminar to turbulent flow. At the very first appearajice of 
amplified vortices the flow still is entirely regular; transition to tur- 
bulent flow may not be expected until the Reynolds numbers are higher. 


1. INTRODUCTION 


Until now the stability calculations of laminar two-dimensional fluid 
flows on straight walls had usually been based upon disturbances in the 
shape of plane wave motions which travel in the direction of the flow. 
After some initial failures (see Noether's comprehensive report, 1921 
(ref. 2)), the reseeirches by Prandtl, Tietjens, Tollmien, and Schlichting 


*"Uber eine dreidimensionale Instabilitat laminarer Grenzschichten 
an konkaven Wanden." Ges. d. Wiss. Gottingen, Nachr. a. d. Math., Bd. 2, 
Nr. 1, 1940 . 

^Presents part of a thesis by H. Gortler submitted in partial ful- 
fillment of the requirements for the degree of Dr. Phil, in Math - Natural 
Sciences Faculty of the University of Gottingen. 
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have, since 1921, produced results which compared well with observations 
and to a certain extent yielded information about the important question 
of the origin of turbulence from small disturbances. Schlichting, 1954 
(ref. 7)> gave a report on the results of these investigations. A brief 
glance at the method is indicated. 

In these calculations, a velocity distribution U(y) that depends 
only on the coordinate y at right angles to the plane of the wall is 
assumed as the basic flow. The omission of variations of the laminar 
basic flow in the x-direction (= principal flow direction parallel to 
wall) was dictated by mathematical reasons; and the results of the calcu- 
lations enabled valuable deductions to be made as long as the variations 
in x-direction were not excessive. To the basic flow U(y) were added 
distxurbances of assumedly sufficient smallness to permit linearization of 
the hydrodynamic equations with regard to the components of the disturb- 
ance. This way the problem could be narrowed down to an expression for 
the stream function of the disturbance in the form 


\lf(x,y,t) = cp(y)e^(°^"P't) 


( 1 . 1 ) 


A particular disturbance can then be built up by the Fourier method as a 
disturbance of a general kind by a linear combination of such partial 
oscillations, While a is assumed as real, the prefix of the imaginary 
part of ^ determines whether there is amplification or damping with 
increasing time t. 

The more general expression of three-dimensional disturbances in the 
form of traveling waves, which are parallel to the flat wall but oblique 
to the base flow direction, hence, for which the velocity compo- 
nents Uj^(i = 1,'2,5) are given by 


^i = fj_(y)ei(a.ix+a2Z-pt) 


( 1 . 2 ) 


(z coordinate parallel to wall and perpendicular to principal flow direc- 
tion), was analyzed by H. B. Squire (ref. 6). By comparison with the 
aforementioned special case (l.l) to be treated Independently, he was able 
to show that, in the case of the disturbances (1.2) with / 0, og 0, 
amplification always occurrs at higher Reynolds numbers than in the case ^ 
of the disturbances (l.l) with ct^ = a;L2 + 0^2. Therefore, the investi- 
gation can be limited to two-dimensional disturbances of the form (l.l). 

The stability investigation of laminar bovuidary layers relative to 
these disturbances was also applied to curved walls (x is then the arc 
length of the wall). Tollmien's claculations for the flat plate with 
allowance for friction were applied by Schlichting (ref. 5) to the case 
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of flow within a rotating circular cylinder. The stabilizing effect of 
the wall curvature is such that the critical Reynolds number, formed with 
the displacement thickness 5* of the boundary layer, increases with 
increasing 67 R (R = radius of circular cylinder). This stabilizing 
effect corresponds likewise with the concepts associated with the action 
of the centrifiigal force (compare Prandtl, ref. k) . 

Boundary-layer flows on slightly ciirved stationary walls were inves- 
tigated by the writer (ref. 10 ) for stability against two-dimensional dis- 
tvu-bances of the form (l.l). Tollmien's result for flat walls, with fric- 
tion neglected and, hence, with the evaliiation of the critical Reynolds 
number disregarded, was the well-known stability criterion which states 
that boundary-layer profiles with inflection point are imstable. Such 
profiles are characterized by a pressure rise from the outside of the 
boundary layer in the direction of flow. For curved walls, this criterion 
is modified to the extent that , instead of the stipulated change of the 
sign of U"(y), a change of sign at U" + ^17' is necessary (curvature 

radius R of wall positive on walls convex to flow, negative on walls 
concave to flow). This Tollmien instability occurs, therefore, on concave 
walls only after the minimum of the pressure impressed on the boundary 
layer from without, and on convex walls already before the pressure mini- 
mum. However, the effect of the wall curvature is extremely small. 

It is surprising that convex stationary walls in this sense act 
amplifying, but concave walls stabilizing, hence, that the effect of the 
centrifugal force does not appear. A confirmation of the criterion fol- 
lows from the fact that the same can be applied also to Schlichting ' s 
case of a rotating cylinder, as explained in detail in the aforementioned 
report. There the criterion, in accord with Schlichting ' s results, yields 
a stabilizing effect of the rotating concave wall. In impublished calcu- 
lations, Schlichting investigated the case of the stationary curved wall 
in analogy to his own and Tollmien's calculations for the flat wall, with 
allowance for friction, for the purpose of observing the wall -curvature 
effect on the critical Reynolds number. In a personal conversation, 

Mr. Schlichting told me that these calculations also proved the stabi- 
lizing effect of concave walls and amplifying effect of convex walls. 

In the present report, it will be shown that boundary-layer profiles 
on concave walls can become unstable relative to certain three-dimensional 
disturbances. It involves an instability that does not occiir on flat or 
convex walls. The friction is duly allowed for in the calculations and 
even the impediment due to the now more complicated type of disturbance 
can be overcome in a relatively simple manner. As to the type of these 
disturbances, they are similar to those investigated by Taylor in 1923 
(ref. 3) in the flow between rotating cylinders and which led to the well- 
known instability (excellently confirmed by experiments Taylor made at the 
same time) in the form of appearance of sharply defined vortices distrib- 
uted boxlike in rectangular zones (compare fig. 1 ) taken from Taylor's 
report . 
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The expression of corresponding equidistant vortices in the boundary- 
layer at a curved wall, ' in which the axes of the vortices are parallel to 
the principal flow direction (see representation in fig. 2 ), leads through 
the Navier-Stokes equations and the equation of continuity to an eigen- 
value problem: The amplification of the vortex disturbance for a pre- 

scribed basic flow on a given wall must be computed for each vortex dis- 
turbance and each Reynolds number of the basic flow; in particular for the 
neutral disturbances (zero amplification) a "critical Reynolds number" for 
each vortex distance must be determined. It is interesting to know how 
these results tie in with the type of basic flow and the wall curvature, 
and also the size of the vortices which are amplified first at increasing 
Reynolds number, as well as the question of the most dangerous vortices 
from the point of view of transition from laminar to turbulent flow. 

At the present state of the experimental investigations, only the 
order of magnitude of the effect is of interest. All the calctilations 
are centered on these claims, hence do not aim at an exhaustive mathemat- 
ical treatment of the present dist-urbance problem but rather to a reply 
to the questions of interest in practice with an expenditure justifiable 
to the claim. 


2. DEVELOPMENT OF THE DIFFERENTIAL EQUATIONS OF DISTURBANCE 


Consider the case of two-dimensional flow of a viscous fluid on a 
slightly curved stationary wall. The finite curvature radius R of the 
wall is, for the sake of simplicity, assumed as constant, and R is 
assumed great compared to the boundary-layer thickness B on the wall 
formed under the influence of the viscosity; R is chosen positive for 
walls concave to the flow - since the instability to be explored occurs 
only on concave walls - and negative for walls convex to the flow. 

The basic flow is along the x-direction (x = arc length along the 
wall), y(i 0) is the vertical wall distance, and' z is the coordinate 
at right angle to both in the direction of the cylinder axis out of whose 
surface portion the wall is formed. 
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In these coordinates, the first Navier -Stokes equation, for exaimple, 
reads in full rigor and generality 



+ R u 


V ^ 

uv , 



St 

-f. - . _ 

R - y 

— + 
dx 

Sy 

R - y 

■u- . 

^z ■ 



.R 1 ^ 

4 - - 

r r 2 b^ . 

£u . 

^2u 

R 

- y P dx 

* V 

Ur- 

y)2 Sx2 

^2 

^z2 


1 ^ _ 2R ^ u 

R - y dy (R - y)2 Sx (R - y)2 


where u, v, w are the velocity components of the total flow in x-, y-, 
and z-directions, p the pressure, p the density, and v the kinematic 
viscosity of the flawing mediimi. All flow variations in x-direction are 
disregarded as customary, and R is assumed great with respect to 6 by 

binomial development of ^ ^ ^ and The Navier-Stokes equa- 

tions and the continuity equation, up to the terms of the order ^ , read 
then ^ 


^+v( + v^=v( _ 1 ^ 

ht \ by B. j bz \ ^ ^ 


^ R Sz P by I ^2 Sz2 R ^ 


!«: + v^+ w^ 
ot ay oz 


i — + v| 1 bw 

P Sz I ^y2 ^2;2 R by 


( 2 . 1 ) 


bv V . ^ ^ 
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The undisturbed flow u = UQ(y,t), v = 0, w = 0, p = Pq, which 

itself is to be a solution of the hydrodynamic equations, for which, 
therefore, 

^Uq _ 1 ^ 

^ R 

"0 1 ^0 

R P Sy 


are applicable, is to change very little during the interval in which the 

disturbances are to be observed. Therefore, and, hence, its equiva- 

ot 

lent viscosity term is deliberately disregarded hereafter and Uq is 
put = UQ(y) . 


This basic flow Uo(y) is a laminar boundary-layer flow formed by 
some previous history based on the viscosity effect. Use is made occa- 
sionally of the conventional idealization of such a boundary layer, which 
consists in assuming instead of the asymptotic transition in the outer 
flow xjq = Uq = const., an increase of Uq(0) = 0 at the wall up to the 

value Uq(6) = Uq at a certain point y = 6 = "boundary-layer thickness," 
while putting Uq = Uq for y ^ &. The minor effect of the assumedly 

slight wall curvature on the outside flow is ignored, since it plays no 
part within the framework o:^ our theory of a first approximation. 

On the assumption that R » 5 

*1 ps. 

the term - — with respect to 
R ay 

tions (2.1) on account of ^ ~ ~ — The same applies to the 

oy & Sy^^ 5 ay 

two other velocity components. The essential effect of the wall curvature 

2 

becomes evident in the term ^ of the second equation (2.1). Moreover, 

no systematic difficulties are encountered if the cited small terms are 
carried along in the subsequent calculation. But, since they only hamper 
the task and contribute nothing to the effect involved, they are 
discounted. 


, the term relative to ^ and 

2 S ¥ 

2-^ can be disregarded in the equa- 
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So, in conformity with the arrangements at the beginning, the fol- 
lowing disturbance equation is used: 


u 

V 

w 

P 


Uo(y) + U2^(y) cos azeP"*^ 
■VQ^Cy) cos aze^^ 

^U^Cy) sin azeP^ 

Po(y) + Px(y) azeP't 


> 


( 2 . 2 ) 


a is to be real and the calculation for p itself is to result in real 

values; a = where A is the wave length of the disturbance. The 
A 

quantity p governs the amplification or damping of the flow, depending 
upon whether it is greater or smaller than zero. The equation (2.2) cor- 
responds to a vortex distribution at the curved wall, the axes of which 
coincide with the direction of the principal flow. Figure 5 represents 
the streamline pattern in a section normal to the principal flow direction. 

Introduction of equation (2.2) in the equations (2.1) following the 
omissions arising from R » 6 results in the linearized equations with 
respect to the disturbance 


^“1 + ''1 ^ H 


d u 


.<iy^ 


1 ^2„ 
— - a u- 


2Uq 1 dpj^ 




“ a‘-v. 


(2.3.1) 

( 2 . 5 . 2 ) 


a 2 

P^i - p Px = ^(^2 


(2.3.3) 


1 ^ 
a dy 


(2.3.i^) 


They apply as long as the disturbance velocities are small with 
respect to the basic-flow velocity. 

To treat this system of ordinary differential equations for the 
unknown functions Ux, Vx, Wx, and px, we insert Wx from (2.3.4) in 
( 2 . 3 . 3 ). The result is p^ as a differential expression of the third 



8 


NACA TM 1575 


order in v^. On substituting this expression for in (2.3.2), u^ 

appears as differential expression of the fourth order in v^. Combined 

with (2.3.1), the following system of coupled differential equations is 
obtained for u^ and v^: 


V 


d^m 


(p + va^) 



(2.i^.l) 



d^v 


(p + 2va^) + a^(p + va^)v^ 


2a^UQ 

u, 

R 1 


(2.4.2) 


When Uq^ and vq 

and (2.3.5). 


are known. 


wq and p^ are computed from (2.5.4) 


It is not recommended to set up a differential equation of the sixth 
order for u^^ or for alone by further elimination. The subsequent 

calculations are rather based direct on the systems (2.4.1) and (2.4.2) and 
merely produce a simplified mode of writing. With 5 denoting a suitably 
chosen measure for the boundary-layer thickness, the following dimension- 
less factors are utilized: 



(2.5) 


For neutral distvirbances, that is, that state of transition in which 
the distxirbances are neither amplified nor damped, P = 0, hence t = a. 
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It further is appropriate to use the quantities 


u' = 


Uo&\-i 


u. 


v' = Vn 


( 2 . 6 ) 


instead of u-j_ and The prime is also omitted in the following with- 

out running a chance of causing a mixup with u and v defined by (2.2). 


The differential equations (2.4.1) and (2.4.2) can be written briefly 
as differential equations for u and v, as follows: 




, dU 
Lu = — V 
dTi 

LqLv = - o^iiUu 




(2.7) 


by utilizing the differential operators 


L = 



L 


0 



In conformity with the order of this system, six boundary conditions 
can be prescribed. It is especially stipulated that ui(0) = vi(0) = 

■W]_(0) = i.e., that the fluid hugs the wall. So with consideration to 

(2.3.4), it is required that u(0) = v(0) = v'(0) = 0. With the other 
three conditions, the decay of the disturbance at t| — > m is attainable, 
or when the boundary layer at y = 6, that is, t) = 1 is permitted to 
change to the cons'tant ou^side flow. Thus the smooth junction of three 
disturbance components with the respective values decaying with q — > <x> 
outside the boundary layer is assured. The symbol "<»" signifies 
"sufficiently great." 

The homogeneous system of differential equations (2.7) together with 
six homogeneous boundary conditions produce an eignevalue problem for the 
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proposed values of U(t]) and R/ 6 : The magnitude of amplification p 

for every given wave length A and every given Reynolds number Re = 

must be determined (i.e., the relationship existing between the parameters 
T, o, and n, required for solving the homogeneous boundary value prob- 
lem, must be calculated). The neutral disturbances (3 = 0, that is, 

T = a) especially, call for the determination of a "critical" Reynolds 
number of every wave length of disturbance A, at which the particular 
disturbance is exactly maintained without amplifying or decaying. 


In the subsequent analysis of the eigenvalue, the practical aspect 
is the primary object — namely, at what Reynolds number does amplification 
appear (stability limit)? What is the range of the wave lengths of dis- 
turbances that can be amplified at all? At what wave lengths does ampli- 
fication appear first when Re increases? What disturbances are ampli- 
fied most and are therefore most dangerous from the point of view of 
turbulence? What effect has the amount of the wall curvature on these 
data? Are there appreciable differences when different boundary-layer 
profiles U(t|) are used as basis? The question of calculating the eigen- 
function is disregarded in the present report, although it may be stated 
that the method developed enables an approximate representation of it. 


It is readily apparent from (2.7) that the Reynolds number and the 


wall cvirvature appear only in the form of the dimensionless 
(namely, in parameter n). 


V 


5. CONVERSION OF THE DIFFERENTIAL EQUATIONS OF DISTURBANCE 
TO AN EQUIVALENT SYSTEM OF INTEGRAL EQUATIONS 

is identified by the following postiilates 
at T) / is twice differentiable with 

T] tIq in 0 ^ ^ CO. 

at the point q =11^, but has in its 
first derivative the discontinuity defined by 

0 t (’'0 • ’'0)} ' - ^ 

( 4 ) G(^; = 0 and G also disappears at t) — 5»<». 


Green's function G^q; 

(1) G^q; qQ^ in 0 g q < 00 

respect to q . 

(2) LG s ^ G = 0 at 

dq2 

(5) G^q; qQ\ is continuous 
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Green's function tiq^ is to have the following quality: 

(l) is four times continuously differentiable with 

respect to t) in 0 ^ H ^ » at y\ -f tiq. 


(2) LqLH = 0 at T ilO 0 ^ rj ^ co. 

(3) At the point ti = t)q, is continuous including first 

and second derivatives, but the third derivative has the 
discontinuity 


llm 
e — ^0 


d^H 


f (no + e; no) - H (^0 - no) 


dq 


dq' 


= - 1 


(4) H(0; tIq) = ^ (O; t^q) = 0 and H disappears at q — ^ 

^ these requirements, G and H at 0 < q < oo are clearly 
identified. The calculation gives 


G(n; no) = 


(a+T)(a-T) ^ 


i e“''’no sinh rq for q ^ q^ 
e“''’n sinh rqQ for qQ ^ q 

; _ e'^^O^^cosh oq - cosh rqj 


H(n; no) = 


- ^Te~^0 - ae sinh 


oq - a sinh rq) 1 for q ^ q 


0 


(a + t)(o - t)' 


^ e"*^ - e“^n^ (cosh 


«no 


lTe“°^ - ffe"^^ H T sinh oqQ - a sinh 


CT 


- cosh 


Tno 


^no) 


for qQ ^ q 


(3.1) 
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In the event that 3 = 0, becomes 



The differential equation system (2.7) is equivalent to the integral 
equation system 

uM . g(„; „o) 

* (5.2) 

v(q) = o2n H(n; 

4. METHOD OF DEFINING THE EIGENVALUES 

To begin with, the integration interval is divided into partial 
Intervals of the same length d, and signify points 

of the k-th partial interval: (k - 1) d ^ kd (k = 1, 2, 3 .-.), 

and where, for simplicity sake, q (k) ^ 

The subscript k added 

to a function symbol Indicates that the particular function is to be 
formed at a point of the k-th partial interval, say about 

furthermore, Gj^j^ = G^q^^); %k = ^ reason of 

the symmetry of the Green functions, G^j^ = and H^it = Hjj.j^. 
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Patterned after the Fredholm theory the integral in (3.2) is 
replaced by summation 


«i 


dZGikUk'vk 

k=l 


(4.1a) 


Vi = o^^xd^HikUi^Uj^ 

k=l 


(i = 1,2,5,...) 

(4.1b) 


Ujj^' disappears for sufficiently great arguments ilQ(k). Letting, as 

approximation to the asymptotic transition, the boundary layer change to 
U = Constant at H = 1 (that is, y = 5), it can be stated more accu- 
rately that = 0 ^ ^us at finite d a finite sum is 

involved in the summation (4.1a). If nd = 1, the summation along k 
must be extended from 1 to n. As a result, only the with k ^ n 

appear on the right-hand side of the equation (4.1a). 

Correspondingly, considering only the equations with the Vj^ at 

which i n in (4.1b), the infinite sum on the right-hand side can be 
approximately replaced by a finite simi of k - 1 up to a sufficiently 
great k = N, because the values decrease rapidly with increasing 

k owing to the upwardly restricted 1 ^ n (the Uj^ themselves decay 

with increasing k). The homogeneous system of the n + N equations is 
therefore investigated 

k=l 
N 

o^qd^I^HikUkUk - Vi = 0 
k=l 

for the N unknown Ui(i = 1,2,...N) 

Vi(i = 1,2, . . .n) . 


(i = 1,2,. ..,N) 

(i = 1,2, ...,n) 


r (^- 2 ) 


and the n unknown 
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The vanishing of the determinants 


1 

0 

0 

GiiU^ 

d 

Gi^U^-d 

•“ ^InV 

0 

1 

0 

G U ' 
21 1 

d 

G U 'd 
22 2 

... G U ' 

2n n 

0 

0 

1 

^nA 

d 

^2^2’*^ 

... G-, U ’ 

Nn n 

o^pHj^lUid 

o^pHq^^U^^ ... 


- 1 


0 

0 

o2pH22^U^d 

aSpHggUgd . . . 


0 


- 1 

... 0 

a2pHniUid 

c2ni^U2d . . . 

u2pH^Ujjd 

» • • 

0 


0 

• ••••• 

- 1 




postulated for the existence of a nontrivial solution system leaves an 
algebraic equation for p at given U, a, and t and, especially, the 

critical value of the dimensionless for t = a. 

V \1R 


Every point within the boundairy layer contributes two non- 

trivial series to the determinant — that is, two series (or gaps) in which 
not only the terms of the principal diagonal are different from zero; a 
point at the border or outside of the boxmdary layer supplies only a non- 
trivial series; the wall point (say, chosen as produces only 

trivial series because there the Green functions (and U also) disappeair. 


5. CALCULATION OF CRITICAL REYNOLDS NUMBER VARIATION 
IN ITS MINIMUM WAVE LENGTH RANGE 


A few words concerning the choice of basic flow for the proposed 
nvunerical calculations are indicated. Theoretically, the basic flow U 
represents any boundary-layer flow formed at a wall due to friction and 
some earlier history. The present calculations are based on the data of 
the Blasius boundary layer of the flat plate (ref. 8). Several other 
profile forms are included for comparison. 

As regards the profile U of the plate boundary layer, the wall 
distance at which the boundary layer in its asymptotic transition to the 
outside flow diverges only 1 percent from this flow (curve 1, fig. 5) 
serves as measure 6 for the boundary- layer thickness. This is the wall 

distance at which the variable 
value 5. 


—\l— in Blasius 's report assumes the 
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As a practical check on the quality of convergence of the calcula- 
tion method developed in section 4, ti( 4) _ n^(k) = ^ ~ ^ and d = i 

' ‘On n 

were selected and the following three approximations calciilated; for 
q(k) and points 0, l/2, 1; 0, l/5, 2/3, 1; 0, l/4, l/2, 

5/4, 1 were taken. According to the remarks made at the conclusion of 
the preceding section 4, the calculation of three-, five-, or seven-row 
determinants is involved, which result in a linear or quadratic or cubic 
equation for with respect to o and t. Evaluation for neutral 
disturbances (t = ct) showed that, to each value of the parameter a = oS, 
that is, to each wave length of disturbance there corresponds the related 
value of n as smallest root. Since the equations exhibit, on the whole, 
coefficients with alternating prefix, only positive roots p. are obtained 
by this calculation, that is, positive values of the critical dimension- 
less p, hence an instability of the assiomed type only on concave 

walls (R > 0) result. 

The results of this preliminary calculation are shown in figure 4. 

The convergence for the parts of the curve above greater or smaller oB 
values, where the curves continue to rise, was not quite satisfactory. 

But the range of the minimum, which is of chief interest here, emerges 
sufficiently accurate. 

Beyond these approximations, other points T](k) outside the boundary 
layer were assumed for individual oB values as a check that the 
approximations achieved in figure 4 are not subjected to appreciable 
changes. The minimum becomes a few percent less and shifts slightly 
toward smaller oB values. 

Incidentally, it should be noted that the order of magnitude of 
these numerical values had been checked by special calculations. Origi- 
nally it had been attempted to solve (2.7) by expanding t\ in power 
series. The convergence for u and v from t] = 0 on was very slow. 
Therefore, series from = 1 on were resorted to. Corresponding to 
the order of the differential-equation system and the number of boundary 
conditions, three coefficients each had to be determined. In consequence, 
u and V had to be joined continuously with continuous first and second 
derivatives within the boundary layer. This gave six linear homogeneous 
equations for the six still indeterminate coefficients and the stipulated 
disappearance of the six-row determinant of this equation system produced 
the conditional equation between n, a, and t. But these calculations 
failed at the evaluation of the determinants. The values of u pnd v to 
be gained from the series and their derivatives could still be determined 
with an accuracy of 1 percent at q = 0.5, but on account of the lanavoid- 
able large figures appearing in the solution of the determinants, the 
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results could no longer be regarded as reliable. On the other hand, 

Un5 

near the minimum on the curve of against o6, they yielded 

results which in order of magnitude agreed with the previous calculations 
Because of the surprisingly small values of the critical Reynolds number, 
the new calculations explained above were carried out. 


The next step was to find 
by a different choice of basic 

with the T|(^) places 0, l/4. 


layer profiles 




U(ti) 


Sin 


n n + e 

2 1 + G 


the extent of the change in the results 
flow U(t]). To this end the calculations 

1/2, 3/4, 1 were repeated for the boundary 


sin 


€ 

2 1 + e 


)( 


1 


sin 


rt e 
2 1 + e 





for 0 T| ^ 1 
1 for T) > 1 


r (5.1) 


with ^ ^ ^ = tl, that is, = 1.752, ^2 ^ ~ 0-5890. The first pro- 

file (e = 1 ) has negative curvature throughout, the second, (e = Cg) has 


an inversion point. (See fig. 5, curves (2) and (5).) 


To assure a physically logical comparison of the results for the 
several boundary-layer profiles, it was postulated that all profiles have 
the same momentum thickness 


•a 



Uo)uQdy 



- U)UdTi 


( 5 . 2 ) 


which is a measure for the loss of momentum in the boundary layer. This 
condition is met when between the individual boundary-layer thicknesses 
the relation a = 0.1116 = 0-13262^ = O.I 3782 exists. Here 6 denotes 

the previously defined thickness of the Blasius plate boundary layer, 

61 and 62 the boundary-layer thickness (ti = 1 ) for the sine profiles 
( 5 . 1 ) with ex and €2. 

The resiilt of the comparison is shown in figure 6. It was found 

u a 

that, when ~^\l^ plotted against cl‘8, the individual cirrves within 
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the scope of our approximation do not differ appreciably from one another. 
(a corresponding comparison based on the displacement thickness 6* 
instead of -a produces curves which differ from one another considerably.) 


A final calculation, as the roughest approximation to an actual 
boundary-layer profile, was made on the section profile 


U = 


jrj for 0 ”2 T) < l| 
|l for T) ^ 1 


(5.5) 


(compaxe curve 4, fig. 5)* If 


6^ is the boundary-layer thickness of 


At identical momentum 


this profile (y = d-z at tj = 1), then d = ^ 6- 

thickness with that of the profiles used so far, the difference is 
slightly greater, but, considering the rough approximation (5.3), the 
departure from the results so fair is not very great. (Compare curve 4 
fig. 6.) ^ 


The amplifications in the explored wave-length range, at least in 
vicinity of the critical Reynolds numbers, can be determined by' the same 
approximte method. These calculations were made on the Blasius plate 
boundary-layer profile. Instead of the extreme case (3.1a) of Green's 
function H(t]; Tq) ; 4he more general expression is obtained from (3.1). 

362 

that is, to each pair of parameters a, t, there 


To each oS and 


corresponds a particular value of the dimensionless 

V 


V 


— . The curves 
R 


= Constant are obtained by graphical interpolation after conversion 

2 

of 5 to i3. (Compare fig. 8.) For greater parameter values the 

quality of the approximate calculation decreases quickly. 


6. ASYMPTOTIC STATEMENTS 


Supplemental to these results for great and small values of oS, 
a few statements are indicated. A differential equation of the sixth 
order for u alone can be obtained from (2.7) by elimination of v. Its 
form is disagreeable for the general calculation, but it enables a pre- 
diction to be made for the extreme cases of great and small values of a 
and T. In this differential equation the coefficients relative to a 
and T represent polynomials up to the sixth degree. Considering only 
the two highest powers on the assumption of sufficiently great values of 
a and t inside the boundary layer, the problem reduces to the second 
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order differential equation 

r2(2a2 + t2)U'2u" - 2 t 2 ( ct 2 + t2)u'U"u' + 

+ t 2 )( 2 U "2 - U'U'") - a 2 r^U' 2 ^ u = - |iUU' 5 u ( 6 . 1 ) 

Integration of this equation across the boundary layer gives the 
relation 


Uq6 2 5 


V R 


J - T (o^ + t^)( 2U"^ - U'U''')^u + 2 t^(o^ + T^)U'U"u' - t^(2o^ + t^)U'^u" 

2cj 2 r UU'^udTi 

^ n 


dti 


( 6.2 


The integrals still contain the unknown function u and its first and 
second derivative, but the derivatives multiplied by polynomials of lower 
degree in a, t, so that the essential contributions to the integral are 
already included in the estimation 


V\26 

. V / R 


2 



U'^udT] 


UU'^udTi 


(6.2a) 


Equation (6.2a) is evaluated by an approximation expression for u by 
means of a polynomial of the fourth degree in r\, taking into considera- 
tion the boundary conditions u(0) = 0, u"(0) = 0 (hence v(0) = 0), 

u"'(0) - t 2 u'( 0 ) = 0 (hence w(0) = 0) and u' (l) + tu( 1 ) = 0 (constant 
connection with constant tangent to the solution for 'u outside the bound- 
ary layer, which according to (2.7) is given by u = Constant e-TT] on 

account of U'(ti) =0 for t] ^ 1 and hence u" - t 2 u = 0 ). As a result 
u is closely approximated in wall proximity and the postulated decay 
toward the outside is attained. Minor errors in u near the outer edge 
of the boundary layer are of no consequence in view of the rapidly 
decreasing U'; errors of u in the numerator and denominator act in the 
same direction (positive integrands throughout), thus affecting the result 
very little. Again the asymptotic relation (6.2a) manifests the existence 
of instability at concave walls only (R > 0 ) . 


The defined polynomial for u reads rigorously 


u = constant + 1 ^)ti + t^{t + 4 )^ _ ^1 + t + ^ + 
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but in practice only the highest powers of t are effective for great 
T. Figure 7 represents this approximate function u for several values 
of T . 


The evaluation of the above appraisal for very great cr and x 
gives for the Blasius plate profile the asymptotic formula 

(6.M 

corrected for d . , 

The same calculation for the section profile (5.5) gives the factor 
2.1 instead of 2.5, hence^ a slight difference only. The tie-in with the 
results obtained for average a values is readily accomplished with the 

asymptotic formula (6.4). Figure 8 represents the variation of the crit- 
Uqi9 rr 

ical factor — _ plotted against a-6 in the double logarithmic net 

2 2 
(curve = 0). The first amplification curves = c = Constant > 0 

are also shown. The variation of these curves at high oi9 values is 

Uqi9 fY 

obtained by addition of 2.5c to the critical values of ~^\/p equal 
o6, as is readily apparent from (6.4). Moreover, by (6.4) 

^ ^ - 0 . 43 ^ ( 6 . 5 ) 


at great cvS, which constitutes an upper limit for the dimensionless 

amplification quantity M solely dependent on ii-. 

Uq R 

An asymptotic prediction for small a is obtained also by an appro- 
priate analysis. It is found that the critical factor increases 

proportional to (od) ^ with decreasing od, as expressed in figure 8. 
(For a more accurate prediction, data about the sixth derivative of u 
and V are necessary.) 


7. DISCUSSION OF THE RESULTS 


On the basis of the data collected in the foregoing, the questions 
formulated above can now be answered in some detail. As regards the 
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stability limit, that is, the Reynolds nijmber at which vortices of the 
particular type can exist for the first time without decaying again, is 


M 

V 


0.58 





(7.1) 


It involves vortices at which ad « 0.l4, that is, whose wave length 
A is given by 


A « 45d(= 5. 06) 


(7.2) 


For the Taylor vortices between stationary outside and rotating 
inside cylinder (ref. 5), the vortex appearing at the stability limit has 
a wave length of about double the distance of the two cylinders, the vor- 
tices thus filling quadratic cells (fig. 1). In the present case they 
fill cells with a width of about 2-1/2 times the boundary-layer thickness; 
they even extend beyond the boundary layer. The stability limit for the 

Taylor vortices is given by = 4l.5^ where d is the cylinder 

spacing and Uq the velocity of the rotating inside cylinder while the 

outside cylinder is at rest. At d = — d, = 2. 

6 V 



The appearance of the first vortices in the boundary layer does in 
no way indicate incipient turbulence of the flow. On the contrary, it 
should be emphasized that the flow will be regular in every way, just the 
same as before. (Naturally, it does not include the case in which ordi- 


nary plate turbulence already occurs at very great 



No incipient 


turbulence can be produced until the Reynolds nmbers become considerably 
higher so that the disturbances of an entire range of wave lengths expe- 
rience sufficient amplification. The same holds true for the Taylor vor- 
tices between fixed outside and rotating inside cylinder; the vortices 
first appear as predicted, but the flow does not become turbulent until 
the velocities are higher. 


The theory developed in the present report postulates that the vari- ' 
ation of the flow in principal flow direction is small enough to be disre- 
garded. When the variation in x-direction is small, the results obtained 
retain their validity as good approximations. In consequence it is justi- 
fied, under this hypothesis, to inquire into the fate of a vortex of given 
wave length in its wandering in flow direction through a boundary-layer 
thickening up at constant outside velocity. The momentary shape of the 
boundary layer has no appreciable effect on the results, as already seen, 
when it is referred to the momentum thickness as characteristic length. 
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In the diagram the vortices of constant wave length 

A describe, by virtue of the identity 

(a,)-5/2M^5(a.)5/2 (,.3) 

curves of the configuration 


with (2n)3/2 c = 


Uo>^ 

V 


— = Constant. 
K 


These curves cross the system of 


ciirves of constant amplification and are reproduced for some values of 
UqA rr 

the parameter ^ figixre 8. There are curves in this series which 

cross the zone of unstable disturbances - when they enter it they always 
cross it since the ciirves of constant amplification, and especially the 

curve = 0 at great a values, vary proportional to (cL-d)^ (see 


equation (6.4)) - and there are curves in the series that never reach the 
instability range. Thus, the vortices corresponding to the latter are 
never amplified but always swallowed by the viscosity effect. In the 
extreme case there is a curve which is tangent to the neutral curve 

= 0. This is the case for the carve with the parameter ~^\f^ ~ 5® 

shown in figure 8. Therefore, if the disturbance of the wave length 
in wandering through the thickening boundary layer ever is to reach an 

undamped state, — \/J must be ^ 50 ; that is, the inequality 

V f rv 


- A ^ 13.6/m) 
R ■ V V / 


(7.5) 


must be fulfilled, which affords a measure for the smallest vortices which 
are able to experience amplification at all. At the instant where it 
reaches its solitary neutral state, the particular boundary disturbance 
has a specific wave length referred to the momentum thickness prevailing 
at that point. According to figure 8, the contact of the aforementioned 
cijrves occurs at about a-d w 1.1, where, therefore. 


A « 5.7^3 = 0.636 


(7.6) 
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Therefore, the wave length A of the disturbance must have a certain 
magnitude characterized by ( 7 » 5 ) if it ever is to get in a critical situ- 
ation with the increasing boundary-layer thickness. If equality exists 
in ( 7 - 5 )> this instant is given by the fulfillment of ( 7 * 6 ); damping 
occurs before and after. If inequality exists in (7*5)^ then the critical 
ratio of wave length to boundary-layer thickness is already reached at a 

certain stage, where as yet ^ > 0 . 65 , after which the disturbance is 

amplified until a certain second ratio ^ < 0.63 is reached: from then 

o 

on the disturbance is damped again. However, the prediction about the 
second critical ratio is applicable only when the disturbance on the pre- 
viously transversed path of amplification does not exceed the theoreti- 
cally specified range of "small" disturbances. 


The last question to be answered concerns the most dangerous disturb- 
ances, that is, disturbances in the whole range of wave lengths which in 

traveling through the boundary layer at equal Reynolds number expe- 

rience the highest amplification, or in figvire 8 , the curves 
~^\j^ = Constant, which prevail at the start of their amplification path 
before transition to turbulent flow in the range of minimums of the curves 
— = Constant. 

V 


In an article by M. and F. Clauser (ref. 9), the appearance of tur- 
bulent flow at the concave wall was observed for Re„ = = 2.6 x 10 ^ 

A y 

at point ^ = 0.75 and for Re^ = 5.1 x 10 ^ at point — = O.J+5. Using 
R R 

the Blasius law of growth of the boundary layer at the flat plate 

as basis, the values of the critical dimensionless factor 




^ = 2. m 
5\/uo> 


M 


^1^ are 10.6 and 8 . 6 . They are indicated by the dashed markings in 


figxrre 8 . A rough extrapolation indicates that, in the vicinity of these 
values, the amplification curves = Constant (or the curves 

= Constant = ^ which, based upon the law d = , is the sam^ 

Uq ^ V 3 yUo / 

are minimum at about a = 0.6, hence o5 = 5 . 5 . Assuming that the turbu- 
lence is caused by vortices of the type investigated here, at reversal the 
boundary-layer thickness would, roughly speaking, have increased up to the 
order of magnitude of the width of the highest amplified vortices. The 

latter, in turn, would have a wave length A of about A « 50R^^^ • 
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For a more accurate prediction about the most dangerous vortices, the par- 
ticular total amplification throughout the unstable range would have to 
be determined for different vortices a = Constant by integration; but 

for this, the few amplification curves, which at greater values of 

V 

become unreliable, is insufficient. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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Colored fluid 



■Figure- 1 .- Vortex between the walls of two eoncentric rotating-- 
cylinders according to G. I. Taylor, streamline pattern following 
incipient instability (inside and outside cylinder rotate in same 
direction). 
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Figure 2.- Vortex disturbances in the flow of a fMd on a concave wall, 
axes of vortices parallel to principal flow direction. 
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y 



Figure 3.- Scheme of streamline pattern in a section at right angle to 

the principal flow direction. 
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Figure 4.- The critical factor for Blasius’s flat plate boxmd- 

ary layer plotted against a& computed by three increasing 
approximations. 
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Figure 5.- The boundary-layer profiles of equal momentum thick- 
ness ^ used as basis of the calculation. 
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Figure 7.- Approximate function for u (according to equation (6.3)). 
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Figure 8.- Total variation of the critical factor plotted 

2i 

against aa and the first amplification curves = const. The 

parallel lines in this diagram represent individual vortices 
(A = const.) in a slowly thickened boundary layer at constant out- 
side flow velocity Uq. 
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